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Abstract: In this paper we derive a generic expression, which is valid for scales 
larger than Hubble radius and contains only the local terms, for the second order 
curvature perturbations for more than one held, provided the expansion is sourced 
by the energy density of a single held. As an application, motivated by our previous 
paper [1], we apply our formalism to two helds during preheating, where the inhaton 
oscillations are sourced by potential which is governing the expansion of the 
Universe. A second held a, coupled to the inhaton through is excited from 

the vacuum huctuations. The excited modes of a amplify the super-Hubble isocur¬ 
vature perturbations, which seed the second order curvature perturbations to give 
rise to a signihcantly large non-Gaussianity. Our results show that within 3 inhaton 
oscillations for a range of parameters, 1 < g'^/X < 3, the non-Gaussianity param¬ 
eter becomes: /nl > 0(1000), which is already ruled out by the current WMAP 
observation. 
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1. Introduction 

Inflation stretches the fluctuations outside the Hubble radius which later on enter 
inside the Hubble radius to seed the structure formation [2]. A single held inhation¬ 
ary model generates Gaussian huctuations which is in an excellent agreement with 
the current observations [3]. These Gaussian huctuations are due to the adiabatic 
spectrum of the curvature perturbations [4], for a review see [5], which in a single 
held case is conserved after the wavelength of a perturbed mode becomes larger than 
the Hubble radius For multiple helds there is an iso-curvature component which 

^Perturbations whose wavelength is larger than the Hubble radius are often abused in the liter¬ 
ature as “super-horizon perturbations”, instead they should be called as “super-Hubble perturba¬ 
tions. Note that due to inflation the perturbations are always well inside the causal horizon. We 
particularly thank Robert Brandenberger for stressing this point time and again. 
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can source the curvature perturbations, i.e., as a result the curvature perturbations 
are no longer conserved on large scales. 

Although negligible but a hnite calculable departure from Gaussianity can be 
obtained within a single held inhationary model, with a self interacting inhaton 
potential. However the measure for non-Gaussianity, dehned as /atl (see the precise 
dehnition below), is bounded by the slow roll parameters [6,7], i.e., /nl lh-3e|, 
where the slow roll parameters, e, r/ 1. This trend remains true even for multi 
held case, see for example [8-10]. Unless the slow roll conditions are violated during 
inhation, it is hard to imagine a reasonably large non-Gaussianity can be generated, 
which would really challenge the common lore that slow roll inhation produces mainly 
Gaussian huctuations. For a review on non-gaussianity, see [11]. 

The only exceptional moment is right at the time when inhation ends, when the 
slow roll parameters become large, i.e., e, ?7 ~ 0{1), nevertheless this still leaves 
Jnl ~ C>(1) [12]. Large fjsiL 3> 1, however, can be obtained during preheating after 
inhation as we illustrated in our previous paper [1] 

Particularly a scalar preheating [16] is a phase where after inhation the coherent 
oscillations of the inhaton, tp, dominates the energy density of the Universe and 
another scalar held, a, coupled to the inhaton through, is excited from the 

vacuum huctuations. In this respect preheating is a non-perturbative phenomena. 
Although preheating does not lead to a complete thermalization [17,18], nevertheless, 
in many models preheating can be just a possibility before the perturbative decay 
of the inhaton, for example see the last reference of [16], if there exists favorable 
conditions such as signihcantly large coupling. In reality it is hard to predict exactly 
how the inhaton couples to the Standard Model degrees of freedom, due to lack of a 
model where the inhaton is not an absolute gauge singlet [19] therefore, it becomes 
difficult to pin down thermalization time scales, etc., which can be tested through 
present cosmological experiments. Nevertheless, very recently as we pointed out that 
preheating can give rise to a large non-Gaussianity which can lead to constraining 
certain parameter regions, such as a coupling constant [1,13,14], we regard this as a 
window of opportunity to test the initial stage of preheating. 

In order to understand why preheating is the test bed for a large primordial 
non-Gaussianity. Let us hrst note that in general we require: 

• Signihcant interaction which leads to curved potential that breaks slow roll 
conditions. 

• Some small amplitude isocurvature huctuations which could be amplihed later 
on via non-adiabatic evolution. 

^See also other examples considered in Refs. [13,14] and [15]. 

^An exceptional case of assisted inflation [20] with gauge invariant supersymmetric flat direc¬ 
tions [19] can possibly address this crisis [21]. 



During preheating it is possible to have both the criteria fulhlled. The hrst 
one can be obtained mainly because during preheating the slow roll conditions are 
violated. The second condition is also satished because the a held, which carries 
isocurvature huctuations ^ , undergoes non-adiabatic evolution. The quantum modes 
of a held are excited because of the time dependence on the ehective mass of a held 
(due to coupling between (p and a where p is oscillating) changes sign, therefore 
violates adiabatic evolution of the vacuum. 

Still a reader might wonder why amplifying Sak leads to a large non-Gaussianity? 
After all non-Gaussianity is related to the second order curvature huctuations. Two 
specihc questions could be: 

• Gan a sub-Hubble process inhuence the super-Hubble second order curvature 
perturbations? 

• Gan existing super-Hubble modes amplify themselves to inhuence the second 
order curvature perturbations? 

Note that the second order perturbations take into account of the interactions in 
leading order of the perturbation theory, namely the quadratic combinations of the 
hrst order perturbations act as sources for the second order. In momentum space 
the sources become convolutions and thus couple diherent scales. To answer the hrst 
question, we should also keep in mind that inhation makes the causal horizon expo¬ 
nentially large compared to the Hubble radius. Therefore causality allows various 
modes to mix, i.e., super-Hubble modes in second order can be inhuenced by either 
convoluting the hrst order sub-Hubble modes or the hrst order super-Hubble modes. 
In paper [1] we studied preheating related to m^^p^-potential for the inhaton where 
the resonance ahects the super-Hubble modes in second order through sub-Hubble 
perturbations, i.e., it relates to the hrst question. 

In this paper we provide a generic expression for the second order curvature 
perturbation which are fed by the super-Hubble huctuations in the hrst order. As an 
example we consider preheating in Ay^^-potential for the inhaton where the second 
order super-Hubble modes are ahected. 

Let us now highlight the role of isocurvature huctuations. It is a well known fact 
that the isocurvature huctuations source the adiabatic huctuations and therefore the 
curvature perturbations. The adiabatic perturbations cannot increase by themselves 
on super-Hubble scales without iso-curvature perturbations. This is true in hrst 
order perturbations, in second order case there are additional sources which feed the 
second order metric and curvature perturbations, we shall see them below. Without 

^The fluctuations in a field will be scale invariant and the amplitude of the fluctuations in a will 
be similar to that of the inflaton, i.e., 5a ~ Hinf /2 tt (during inflation and after the first few inflaton 
oscillations). With these comparatively large initial amplitude and scale invariant fluctuations it is 
possible to amplify the fluctuations in 5ak on super Hubble scales. 



having isocurvature fluctuations it is therefore not possible to amplify the first order 
and subsequently the second order curvature perturbations. However note that non- 
Gaussianity is a ratio between the second order curvature perturbation with respect 
to square of the first order curvature perturbations. If both first and second order 
curvature perturbations grow due to large isocurvature fluctuations then the ratio 
need not be large always. 

Therefore in order to get large non-Gaussianity we need a setup where only the 
second order curvature perturbations grow while the first order curvature perturba¬ 
tions do not grow at all. Gan this happen in reality? 

Yes indeed, this can happen when the isocurvature fluctuations do not seed the 
first order perturbations at all, but it only seeds the second order perturbations. We 
achieved this modest goal in our first paper, [1], where the background motion of a 
a held were absent, only the coherent oscillations of ip held were present. In which 
case the perturbations in a decouple from the first order perturbations, but not in 
the second order and in higher order perturbation theory. Therefore exponentially 
large growth in 6ak can feed the second order curvature perturbations to give rise to 
a very large non-Gaussianity. 

The aim of the present paper is to lay down the formalism for second order 
curvature perturbations, where there is a single held dominating the energy density 
of the Universe while other multi-fields can leave their imprints via fluctuations 
We also show that the second order curvature perturbation becomes free from the 
non-local contributions on super-Hubble scales when we neglect the gradient terms, 
which is true in a generic background (all the scalar fields can have non-vanishing 
VEVs), see the detailed discussion in |A.2| . 

As an application we will consider a massless preheating, where the a held does 
not have a bare mass term and the inflaton potential is quartic, i.e., V{ip,a) ~ 
-|- g'^ip‘^(j‘^. This potential has certain nice properties which we will highlight in 
the course of our discussion, i.e., the instability band for the mode depends 
on the ratio, f A, cosmologically interesting solutions can be obtained for the first 

instability band, 1 < g'^/\ < 3. In order to obtain the results analytically, we assume 
that the background value of a is vanishing. A non-vanishing case has to be studied 
separately with an aid of numerics 

A large non-Gaussianity also depends on the number of oscillations, if the number 
of oscillations grow the non-Gaussianity grows exponentially at an alarming rate, see 

®In our previous papers [1,13,14] we mainly studied the first order and second order metric 
perturbations, which need not be conserved on large scales. Here our formalism is more robust 
as we study the perturbations in terms of conserved quantities, Ur the definition see 

below. In this paper we verify our earlier claim that preheating can really boost large primordial 
non-Gaussianity. 

®If a field is also oscillating and contributing to the energy density of the Universe then again 
we would have to address the problem numerically. 
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the Plots. [^ 3 . A simple explanation can be given as follows, flnctnations of the a 
held grow as ~ where /i ~ 0.2 is a physical parameter governing the prodnction 
of a particles ( for a detailed discussion, see below) and x is time in dimensionless 
units which measures the period of oscillations, T ~ 7, so that the time elapsed is the 
number of oscillations, ~ 3, times the period, i.e., x = NT. With these numbers 
the growth factor is amazingly large, 

~ ~ 4000 . (1.1) 

Barring accidental cancellations between various source terms or dampening the ini¬ 
tial amplitude for the isocurvature huctuations this leads to a large non-Gaussianity. 

However we should remind the readers that the phase of preheating cannot last 
forever. Besides the metric and curvature perturbations there is a backreaction 
on energy momentum components due to exponentially large particle creation, this 
inevitably stops preheating. Depending on the couplings the backreaction can stop 
preheating in dozens of oscillations. When the backreaction becomes important our 
assumption with a vanishing VEV of a held also breaks down. 

Unless the time evolution of a held really brings down the non-Gaussianity pa¬ 
rameter to the level that is consistent with the current observations, the model 
considered in our paper with 1 < /\ < 3 is deemed to be ruled out. For other in¬ 
stability bands the ehective mass for the a held becomes heavy, i.e., e// > Hinf, 

which decreases the amplitude of the isocurvature huctuations and leads to a smaller 
non-Gaussianity parameter. In this respect non-Gaussianity acts as a discriminatory 
and it should be considered as a useful observable to distinguish various physical 
situations. 

We begin with the dehnitions, then we describe the hrst and second order curva¬ 
ture perturbations. We then discuss the background evolution where we obtain the 
hnal expression for the second order curvature perturbations for more than one held, 
provided the background evolution is governed by one held alone, for an example this 
could be the inhaton. We then turn our attention to a massless preheating and derive 
an expression for the non-Gaussianity parameter. In order to be self consistent, we 
added an appendix, where useful relationships have been derived explicitly. 


2. Basic equations 

2.1 Definitions 

The metric tensor, including only the scalar perturbations upto the second order, is 
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given by [22] 


( 2 . 1 ) 

( 2 . 2 ) 


9m = (1 + , 


90. = I 


.( 2 ) 


9 .; = o- (1 - 2V-“> - ■!/.«) 5.; + 2B<‘> + 

The energy-momentum tensor for multiple scalar fields is given by: 

1 


T,. = Y. 




0.3 

s9^lv9 '^k,a'^k,f3 




(2.3) 


(2.4) 


All the variables are split into background plus perturbations: 

T = -|- , 


(2.5) 


The evolution for background and perturbations is given by the Einstein equations, 
Gfj,iy = SnGT^u, and the Klein-Gordon equations for the scalar helds, see |A.l . 


Let us specialize to a case when there is a single scalar held determining the 
background evolution, ip, and the rest have constant background values a] =0 
(usually = 0, but this is not necessary). This requires that V^. = 0 through the 
Klein-Gordon equation. The scalar helds are divided as usual into background and 
perturbations: 


(0) 

= ^k 


where is constant and the held huctuations in the hat gauge, = 0, 

are Sasaki-Mukhanov variables, <y 9 d’ 2 ) = and 


+ , 

(2.6) 

+ ^ > 

(2.7) 


2.2 Curvature perturbation 
2.2.1 First order 

The hrst order curvature perturbation is dehned as 


= ( 2 . 8 ) 

where V’ is the perturbation of the trace of the spatial metric and p is the total 
energy density. This combination is gauge invariant. The evolution of the hrst order 
curvature perturbation on large scales is given by, see for example [5]: 


^( 1 )' 


p(0) _j_ p(0) 


(2.9) 
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where is the iso-curvature perturbation or non-adiabatic pressure, 


( 2 . 10 ) 

For adiabatic perturbations, F*^^) = 0, and the curvature perturbation is conserved 
on large scales. This happens when there is a single fluid with a deflnite equation 
of state or for one scalar held on super-Hubble scales. However, in the general case 
where there are multiple fluids or scalar fields, the iso-curvature perturbation does 
not necessarily vanish. 

2.2.2 Second order 

In the second order a similar equation is obeyed on large scales,^ see also [23] 

( 2 )/ ^ n ( 2 ) 

^ p(0) p(0) 


( 2 . 11 ) 


where 




-^(2) - + 2Ff-^ P 


P 


P' 


(0) 


,2 


+ 2 


P 


( 1 ) 


( 0 )' 


P 






P' 


( 0 )' 


v! + 2H^ -n 


p 


( 0 )' 


p 


( 0 )' 


( 2 . 12 ) 


and 


r(2) = p(2) 


„( 0 )' 
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( 0 )' 


p 


( 1 ) 
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( 1 ) 


do)' 


P 


( 0 )' 


p 


( 0 )' 


p 


( 0 )' 


p' 


( 0 )' 


p 


( 1 )' 


p 


( 1 )' 


(0)' \ p(0)' 


p 


( 0 )' 


p(0) -|- p(0) 


F(b% 4 ^(i)F(b _ (2.13) 


Again if isocurvature contribution vanishes then is conserved on super-Hubble 
scales. This has already been shown in [23,25] and it has also been generalized to 
all orders [24,26]. Going to higher orders gives new possibilities, namely could 
be conserved but might not, since it may happen that F^) = 0 but F*^^^ 7 ^ 0 . 
This is especially interesting situation for the cosmological non-Gaussianity which is 
defined as 

= (2.14) 

5 

If Fd) = 0 then (^d) is conserved but F^) 7 ^ 0 then (^d) changes and therefore also 
the non-Gaussianity parameter, /at^, changes. 

^This has been generalized to all orders and scales by Langlois and Vernizzi [24]. 
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2.3 Evolution of the Perturbations 


From now on we adopt the flat gauge, "01 = V ’2 = 0, and neglect terms that contain 
two or more spatial derivatives. We also take only one of the flelds to be evolving at 
the background and call this (p and the rest of the flelds do not evolve at the back¬ 
ground level, cr' = 0. For details and general expressions see [22], The components 
of the Einstein tensor and energy momentum tensor have been given in |A.1| . 

The evolution equations for the background are given by the Friedmann equation, 

= 1 <^(o)'%aVo, (2.15) 

and the Klein-Gordon equation, 

" + 27f(p(°)' + = 0 . (2.16) 


In addition one can write the equation for the evolution of the Hubble parameter, 
which is not independent of the other two equations: 


n' = n‘^- 


2Mj, 




( 0 ) 


/ 2 


The Klein-Gordon equations at 1st order read 

Q™ " + 2HQ«' - v><°> V”' + = 0, 

+ 2HQP' + a^Y. - AQi" = 0. 


Oh 


(2.17) 


(2.18) 

(2.19) 


We will not need the 2nd order Klein-Gordon equations at all. This is due to the 
fact that on the background level only has a non-trivial evolution. 

The gauge-invariant curvature perturbation at 1st order in terms of the Mukhanov 
variables is given by: 

This is a conserved quantity on large scales (with the approximations made ' = 0). 
Using 0 — i Einstein equation one obtains: 



The corresponding conserved variable at 2nd order is given by [22]: 

2 


n( 0 ) ^ 2 ( 0 ) 

a ( 2 ) ^ P ^(2)' , ^ P 

3Uo</2(o)' 


+ -mrQ 


( 0 )/ 




(0) 


VP 


V(0)V^2 


3Uo^2 ^( 0 )' 

,( 0 ) 


+ 




( 0 ) 


/2 


P'- 


SVoH^ 

,/2 


+ 


nVn 


+ 


a 2 p( 0 ) 


-4 


3Uo 


,2 Z_^^k ^ / ow 2 Vn ' 


3Uo p^^y 


k,n 


( 2 . 21 ) 


( 2 . 22 ) 



















Solving and 0*-^^ from 0 — 0 Einstein equations and then taking the divergence 
of the 2nd order 0 — i part of the Einstein equation and then taking the inverse 
Laplacian, we obtain 



(2.23) 


where 


j = _ j_ / nW' ^ -L. q(i) 'nil) I ^(1) 2 r 2 y 

, ^ g^2y^ ^( 0 ) 




/ 2 ' 


P 


(0) 


a? pW" 




/ 2 


E KnQf^ A- E (Qf ’ ■' 

A: ^ k 


(2.24) 


Then the solution to Eq. (p.23|) is given by: 



,( 0 )' 

?r + 



(2.25) 


It would seem that the integral over the source becomes ill-defined at times when 
(^(0)' — which happens during oscillations. During inflation the background held 
is slow-rolling and this problem is circumvented but after inhation when the inhaton 
oscillates before decaying it is a problem. Since and their derivatives are well 
dehned variables with no singularities, the Klein-Gordon equations guarantee this, 
the only singularity in the source integral comes from which gets multiplied 

by after integration. We can extract that singular behavior by expanding all 
the quantities around the time p*, where = 0. Then 


(^(0)' = (p(°)"(pi) (?7 - rji) + 0{r] - rjif -F ..., (2.26) 


where {'qi) ^ 0 during oscillations and the rest of the terms inside the integral 
are also non-vanishing, but regular at p = r/j, so we can replace them in the most 
singular term by constants. Then 




( 0 )' 


n 


drj- 




(^0)' J'P ^ (P -Pi) dp 


{P - Pi? 


const. 


(2.27) 


f2) 

when T] ^ rji so that is completely regular even during oscillations. Next one 
has to define the constant C in Eq. ( |2.25| ) in terms of some initial condition at 
q = po- Here we have to choose a time such that (p^^^'{qQ) ^ 0, since otherwise the 
homogeneous term vanishes and we are not able to define the constant. Then 


)( 2 ) _ 


( 0 )' 


= 




dq 


n 


'Vo 


p 


( 0 )' 


J, 




(2.28) 
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Now we are solving a definite integral over the source term so that there can be 
several points along the integration where = 0, which upon multiplication will 
not be removed, since it is done only at the hnal time. This is not really dangerous if 
we are able to do the integral in closed form but if one needs to evaluate it numerically 
it is a real complication, see |A.3| for a toy example that clarihes the complications 
related to the oscillations. 

Finally, what we are really interested in is the curvature perturbation, and 

f 2 ) 

not the Sasaki-Mukhanov variable However, the curvature perturbation is 

not well dehned during oscillations since it is obtained from the Sasaki-Mukhanov 
variable by dividing with see Eq. (|2.22|) . This problem occurs even in Erst 

order, see Eq. (|2.20|) . The solution is to consider instead ( or (1 +w)C, where 

1 -|- tc = which has been used in Erst order calculations successfully 

and it also works in second order. Using Eqs. (|2.22|,pT^,|2.25|) can be rewritten 


as 


2n 






di]- 




( 0 )' 




/ 2 


E ^ A-1E oi'Tf 




k.n 


p 


hb'n(i) 


(2.29) 


Similar reasoning as we have given for Qip in Eqs. (p.26| , py77|) ) onwards shows that 
" is well deEned. And actually for points where = 0 it is given by the 


second term in Eq. (^2.291) since the integral term vanishes. Another way to derive 
the result for the case of = 0 is to use it in the original Einstein equations, in 
which case one can solve the Mukhanov variable directly at that particular time and 
obtain the result given by the second term in Eq. (|2.29|) . The non-local term can be 
rewritten as (see |A-2|) 


\ k / k 


(2.30) 


The Enal formula for second order curvature perturbation is 


C(2) = C 


n 




n 


A: 




( 0 )'- 


Eq 


/ 2 




k,n 


p 


)(i)'od) 


(2.31) 


Note that the above expression does not contain any non-local term ^ 

®Our expression improves the one obtained in Ref. [22], since we got rid off the non-local con¬ 
tributions by solving them in terms of local contributions, see [A.2|. 
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3. Massless preheating 


As an application to onr formalism, let ns consider a simple potential, 

+(3.1) 

where (p is the inflaton and a held is conpled to the inhaton with a conpling strength g. 
After inhation the inhaton, p, begins coherent oscillations with an initial amplitnde 
~ 0{l)Mp. Dnring this process it amplihes non-thermally. This happens becanse 
a held sees the time varying mass dne to its conpling to the inhaton. In this paper 
we are interested in stndying the hnctnations of the low momentnm modes of the 
a held. This has been stndied in hrst order pertnrbation theory in several papers 
[27-29]. 

For the above setnp the Friedmann eqnation is given by: 

+ ^a^Xp^ , (3.2) 

here we consider the inhaton oscillations is the main sonrce for expanding the Uni¬ 
verse, and the Klein-Gordon eqnation: 

p" + 21-Lp + a^Xp^ = 0 . (3.3) 

Note that for the pnrpose of illnstration we assnme the VEV of cr = 0 In the 
caicniation it is nsefnl to scale ont the expansion by (p = y/a leading to 

= nxf + ^Xx^ , (3.4) 

x" + Ay=* + ^y = 0. (3.5) 


A simple solntion can be obtained for y dnring oscillations while neglecting the 
expansion rate in Eq. O. for details see [27,28]: 


y(p) = yen 



X = \/Xx V , 


(3.6) 


where y is the amplitnde of the oscillations, cn is a Jacobian elliptic fnnetion and 
xq = 2.44 matches the solntion to the slow-roll solntion at the end of inhation. The 
period of oscillations is T = F(l/4)^/y^ 7.4. Then Eq. m can be rewritten as: 

3Mpa‘^ = -y'^ -|- ^Ay^ = ^Ay"^, (3.7) 

®At initial stages cr = 0 is a very good approximation, as the VEV of y decreases it eventually 
sets a motion for a field due to the backreaction of produced a quanta. In Ref. [30] the authors 
analyze non-linear perturbations when the backreaction is important, however they do not estimate 
the non-Gaussianity parameter in their case. In our case we estimate non-Gaussianity from the 
exponential growth of the linear perturbations of cr field. 
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where in the last equality the fact that energy of x is conserved when expansion is 
neglected as guaranteed by Eq. 


The solution of Eq. (p.7|) gives a ^ rj from 
which follows that a" = 0, so that the approximation is consistent. This also results 
into H = I/tj. 

Now we study the gauge invariant fluctuations of cr, which is described by the 
Mukhnaov variable However in this particular case the Mukhanov variable 

is equivalent to the field fluctuations, 5ak, because the background value of a is 
vanishing. Note that although a held carries isocurvature huctuations, it does not 
source the hrst order curvature perturbations due to its vanishing VEV. Therefore 
in our case the huctuations in a cannot be constrained by the current constraints 
from the isocurvature perturbation. However these huctuations do source the second 
order curvature perturbation and therefore are subject to the constraints from non- 
Gaussianity. 

The huctuations of a can be solved from the perturbed Klein-Gordon equation: 

+ 2W<3W' + aV = 0 (3,8) 

It is again useful to make the re-scaling V® change of variable x = 

\fXXf]- Fourier transforming Eq. (|3.8|), we obtain: 


xl-k + 


+ ^cn^{x — Xq, 
A 


,\!,k = 0. 


(3,9) 


where ' = d/dx and Since for observable scales the dimensionless 

momentum is small, for example for length scale of ~ IMpc it is k 10-5°, it is 
important that the instability band contains also the fc = 0 mode, which is actu¬ 
ally true for Eq. (|3.9|) , see [27]. This also justihes taking k = 0 while estimating 
analytically the non-Gaussianity parameter, for a detailed discussion also see |A.4| ^°. 

Gosmologically the most interesting solutions are the ones with 1 < g‘^/\ < 3, 
since the condition for a to be a light held during inhation is that is mass is less 
than the Hubble parameter, more specihcally ( 795 *'°^ < (3/2)i7. Using the Friedmann 
equation this can be written as /\ < 3/16((p*^°V^p)^- During the slow-roll phase 
6 = Ml/2 {V^/Vof = (Mp/(2y2(^(°)))2 < 1, so that during slow-roll the condition 
is obeyed for 1 < g"^/X < 3. This also guarantees that the primordial huctuations 
(generated during inhation) of a held are scale invariant. 

The second order curvature perturbation can be calculated if in Eq. (p.29|) we 
approximate the with its time average (1 -f w)a?p'^^\ where w = 1/3 for 

radiation dominated universe. Then 


e' = eph) + 


+ 


dg 


(1 -FM;)agp(°)(r7o 

n 


■Qa^\V0)Qa\V0) 


'-qq- 


(1 -I- 


09" ^ + 


(1 -I- tc)a2p(°) 

muo 




+ 


p 


( 0 ) 


QW'qW .(3.10) 


lOr 


^The Floquet index defined in Eq. o has a mild scale dependence for cosmologically inter¬ 
esting scales. 
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We simplify the analysis further by disregarding the oscillatory behavior inside the 
integral, since the first two terms are positive definite quantities there is no problem 
with this and the last term in the integral has an extra Hubble parameter in it so that 
it decays away faster than the first two ones and is therefore less significant. Once the 
amplitude of oscillation is known (although now the amplitude grows exponentially) 
the oscillations can be averaged by adding a factor of ~ 1/2 which corresponds to a 
time average of a square of an oscillating function. We also estimate the real solution 
Eqs. ( |A.33yA)^) , see |A.4| , by an exponentially growing function 




(3.11) 


where the 1/2 is due to averaging over the oscillations and /i(s) is given by Eq. (|A.36|) , 
where s = g'^/X. The growth index /i(s) gets values between 0 and /imax ~ 0.2378. 
Since the expansion has been neglected after re-scaling of the held, we take only the 
leading order contributions in terms of the expansion in Eq. ( p.l0|) . After simplifying 
the resulting expression we get 




_L^2„(,)(r-r„) _ _L ^ f j2„(,)(r-r„) 

X Xq \2 / X 


(3.12) 


where r = ^(?7o)/Q^^^(?7o) is the relative magnitude of the a and cp perturbations 

at the end of inhation. If both helds are light one expects r 1. In reality a becomes 
massive during or after the end of slow-roll period, so that there is a slight damping. 
Second the perturbations can have different slow-roll parameters so that r can be 
mildly scale dependent. 


We illustrate the Jnl parameter dehned in Eq. ( 2.14 ) numerically for x = 25 
which corresponds to roughly 3 inhaton oscillations. In Plots |l|,H, we have drawn 
log[/ArL] with respect to x. Note that the /at^ rises exponentially owing to Eq. (|3.12|) . 
This is because the mode, grows exponentially, i.e., see the discussion above 
Eq. (|l|). The same mode sources the second order curvature perturbation, Eq. (|3.12|) , 
but not the first order, see Eq. ( p.20| ), therefore the f^L parameter is alarmingly large. 


This result exceeds the current bound on non-Gaussianity, —58 < fNL ^ 134 set 
by WMAP [3]. Our results show that massless preheating on the parameter range, 
1 < (yf^/A < 3, is already ruled out. We do not expect that our result will be modihed 
dramatically even if the backreaction due to produced a quanta is taken into account. 


Note that the maximal growth index is obtained at g ‘^/A = 15/8 = 1.875 and not 
g‘^/X = 2 as has usually been assumed in the literature, see for instance the recent 
paper [30]. However the difference between the two cases is negligible as is shown is 
Plot 
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Logio [-fNL] 



Figure 1: The evolution of the non-Gaussianity parameter, Jnl, with different values of 
g'^/X = 1.2, 1.5, 1.875 with red, green and blue lines. g^/X = 1.875 is the case with largest 
characteristic exponent. The end point corresponds to 3 oscillations of the inflaton and at 
that point the largest —/nl = 1382. 


Logic [-Lnl] 



Figure 2: The same as previous figure with g"^/X = 1.875, 2, 2.5 with blue, green and red 
lines. The cases with g‘^/X = 1.875 (blue) and 2 (green) prodnce almost identical results, 
which is not surprising since /Ui .875 ~ 0.238 and g ,2 ~ 0.236. 


4. Conclusion 

Following our previous claim [1], we pointed out in this paper that preheating in 
general can produce very large primordial non-Gaussianity. We substantiate our 
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point by providing a formalism which can take into account of the fluctuations for 
more than one held, therefore accounting for the isocurvature huctuations which can 
source the hrst and second order curvature perturbations. The hnal expression for 
the second order curvature perturbation is given by Eq. (|2.31|) . 

We studied a particular case of massless preheating as an example, where the 
second order curvature perturbations are fed by the the isocurvature huctuations but 
not the hrst order. This gives rise to a large non-Gaussianity which already exceeds 
the currently observed limit. The non-Gaussianity parameter grows exponentially 
and the hnal number depends on the number of oscillations before the backreaction 
kicks in to halt preheating completely, see the Plots Based on our result we 

rule out completely massless preheating on the parameter range, 1 < /\ < 3. 

Our results are encouraging and obviously demands thorough study of non- 
Gaussianity generated during preheating. One of the goal will be to setup a numerical 
code where more than one held can participate in the expansion of the Universe. This 
would dehnitely improve the non-Gaussianity bound set in the current paper. 
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A. Appendix 

A.l Setting the background equations and perturbations 

The background Einstein tensor is 


G 


( 0)0 


3 a 


/2 


a? 


Q{0)i _ 


a'^ a" 




a 


The Einstein tensor at hrst order is given by: 


gw; = 4 4^(1), 

^( 1)0 ^ __2 a '^( i ) 


a 


^(i)y ^ 


V2 


4!^0(i) + 2 




(A.l) 

(A.2) 


(A.3) 

(A.4) 

(A.5) 
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and at second order 


g'® + . (A.6) 

= -4 [-0(2) -4- 0(1), (A.7) 

* la a J ^ 

G<^>- , = 1 [-8 f ^,#,(‘),#,(‘l' + - ,#,<') 0 + ^,#,PI' 

^ I \a a J a 

+4^ +2-0(2)(A-8) 

Then the energy-momentnm tensor of the backgronnd is given by: 

(A.9) 
(A-10) 

At hrst order in perturbations 

= -i 0'»''OW' - (A.ii) 

= (A. 12) 

ct 

t(');. = ^ [4(°)'g«' - ¥^(°)''0(1)) -a2V0] <5T, (A.13) 

and at second order 

= -? - V“’V'‘'qw' - 

+4y“>'",).0>"+(jW'T^Oy'AaV2 , (A. 14) 

k 

rA = -0 (i’'">'C-4^'T'“''O0 + 2(3W'(3W + 2 5;oO''(jm^ (a.i5) 

r®' = 1 h<“''<3W' -44 ,i»i><‘)(j(i)' _ yo)'^^(i) 

CL L 

+4y“>'".).0>"+QW'" + 5^Qy'"-o2K2 y., (A.16) 

k 


where 


K) = V{ak = ) (A.17) 

V7 = (A. 18) 

40 = + E ■ (A.19) 

m,n 
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A.2 Huh, getting rid of the non-local terms! 

The non-local term is a complication that every one would love to get rid off. Here we 
work more generally and allow all the scalar helds to have non-vanishing background 
values. Then the Klein-Gordon equations on large scales (neglecting the gradient 
term) are: 




2Wq1'>' + 




aV 


,kl 


Ml 


pi'' 


2 ( 0 )' ( 0 )' ' 
a Vk 

n 


oh = 0 


(A.20) 


Now we multiply the Klein-Gordon equation with djQk, take the divergence of the 
resulting equation, sum over all helds and operate on it by A“^. The result is 






+A-‘E 

k,l 


a Vki 


Ml 


pi' 


2 ( 0 )' ( 0 )'' 

n 


oSohh = 0. 


(A.21) 


Next we re-write the hrst term 




Ql 


I A -1 


- A-^ E. QE C + 2HA-1E. QE C 


+A-^E 


k,l 


, 2 (0) (0) ' 

^2t/, , _ 1 / « Vk G 

a V,ki 

pi 


n 


oSo!‘’)'' = o. 


(A.22) 


In the second term the spatial derivatives and inverse Laplacian can now be solved (as 
it can be done also in the last term) since the potential is symmetric in its derivatives 
and the same is true of the extra mass term 


A-‘ E. (oh'oS)’T - 5 E» (oh')' + 2 WA-> El (oh'oh)" 


'2 Sfc,/ 




pi 


2 (0)' (0)' ^ '■ 
n 


Ohoh = 0. 


(A.23) 


Next rewrite the second term again: 

A-‘E(oh'oA'^ 


-ifEoh'oh) +)5:oh"oh 


+2WA-‘5:(Qi”'o%i5: 


k,l L 


aVi, - 


My 


2 ( 0 )' ( 0 )'' 
«_El_eL 
n 


ohoh = 0 , 

(A.24) 
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and then nse the hrst order Klein-Gordon eqnation to rewrite it as 






5:QW'Qi‘M +2H5 :q«'q 


k 


= 0 . 


(A.25) 


This can be integrated to yield: 


k 


(A.26) 


where /(x) is a fnnction of space only. For one field Eq. (p.21|) shows that / = 0. 
This is not surprising since for one dimensional system the hrst derivative is always 
proportional to the function itself where the proportionality depends only on time 
as long as we work with the leading order of gradient expansion. For multiple helds 
it is enough if the matrix of proportionality is symmetric and depends only on time. 
Special case of this is the one where multiple helds decouple into a system of one held 
models. During slow-roll inhation the potential derivatives are of order the slow-roll 
parameters Then the Klein-Gordon equation at zeroth order of slow-roll is: 




(A.27) 


with the well-known solution (during inhation a = —l/{Hri) with H = const.) 


The derivative of this solution is proportional to the solution itself with the pro¬ 
portionality given by a time-dependent function at the leading order in large scales, 
(actually in this case the time derivative vanishes in leading order 
but in general it is proportional to slow-roll parameters, however the main point is 
that the proportionality factor in leading order does not depend on k). Then the 
non-local term can be solved using the Fourier-transform 


m 

) - e ( 2 ^ 3 ^", 

ra ^ ^ 

f Sk'QZk ■ (k - k')a 


5Z/™(^)(27r)3’^ J 

f d^k'QW ,k . (k - k')a 


m m 

^^Same would be true for the extra mass term in the general case. 














so that the unknown function in Eq. ( A.26 ) actually vanishes /(x) = 0. This 
result is also valid after inflation since one only needs to specify correspondence of 
non-local term to a local term at one time. In a sense this is just a result of the 
fact that after horizon crossing the modes evolve in the same way. Another way 
of getting the same result would be to use slow-roll to neglect the second order 
time derivatives in the Klein-Gordon equation. Then the first order time derivative 
of the Mukhanov variables would be proportional to a time dependent (but space 
independent) symmetric matrix times the Mukhanov variables. Then one can replace 
the time derivative in the non-local term and solve it completely to obtain the same 
result. One should note that both derivations are only valid in large scales since the 
spatial gradient terms were neglected in the derivations. Therefore for small scales 
the scale dependence should be more important so it is not surprising that there it 
may not be possible to remove the spatial derivatives. 


A.3 Time averaging and passing throngh zeros 

Let us study a similar situation with an example concerning elementary functions 
where the integrals can be done explicitly. We take 


f{x) = C cos(a;) -|- sin(a;) = cos(x) 


which is a solution of a differential equation 




dx 


cos^(a;) 


^ + tan(a:)/ = — 
dx cos[x) 


(A.30) 


(A.31) 


The solution Eq. (|A.30|) is similar to Eq. ( p.25|) . It also has a similar property with 
respect to initial values namely cos(a;o) 7 ^ 0 , where Xq is the initial value, has to be 
obeyed. Then 


3 fi^o) - sin(xo) / 3 , . / ^ ^ ^ 

jix) = - - — r - cos(x) -|- sinlx) = cos(x) 

cos(a;o) 


fixo) ^ r 


cos{xo) 


I dx — 5 —, 
xo cos^a:)] 


(A.32) 


The hrst line shows that there are no problems for the solution at any x but the 
second form has the same problem as was seen in Eq. (p.28|) , namely there can 
be divergences within the range of integration. Although, as long as we are able 
to determine the integral function the divergences are of no consequence, but for 
numerical evaluation they are a complication and have to be regulated. Typically 
the regulation is done by taking a time average of the oscillating function inside 
the integral. However, one should be careful under what situations one should do 
this. For example taking a time average in Eq. (|A.32|) , which would mean replacing 
cos^ X —*• 1/2, would result into a linear growth of f{x), whereas the true solution is 
just oscillating. If the source term on the other hand grows exponentially (multiply 
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sinx with e*) there would be inside the integral, too. Now the time average would 
result correctly into exponential growth of f{x) and the behavior of the amplitude 
can be traced even if the phase of oscillation would remain unknown. A good check 
on whether one gets wrong kind of results by time averaging is to use the solution 
at the points when the oscillating function does vanish, since these can be explicitly 
solved from the original equation. 

Multiply Eq. ([A.31| ) with cosx and then take cosx = 0. As a result /(x) = ±1 
at these points which is reproduced by the solution, Eqs. (|A-30| , [Ar3^ , and is actually 
independent of the initial conditions. The same can be applied to Eq. ( p.23|) after 
multiplying by which would correspond to the original Einstein equation and 
then taking = 0. This procedure would exactly give the behavior of at 
these points. The behavior during time of these special points serves as a check for 
the goodness of taking the time averaging. 

A.4 Super Hubble modes 

The region of study we are interested in is the CMB scales. These scales correspond 
to < 10“^° [29]. Therefore for any practical purposes we can take = 0 as an 
approximation. This special case has the beneht that analytical solution exists [31] 

' /?; ^ ; cn\y, 1/v^)) + 1/^2) x 

F(1 + « — 7,1 + /? — 7;2 — 7: cn'^{y, 1/a/2)), 2 nuj < y < (2n + l)a; 
^(^”■+1)^(0!,/d;7 : cn‘^{y.i 1/a/2)) — cn{y, l/\/ 2 )x 

F(1 + a — 7,1 + /3 — 7;2 — 7: cn‘^{y, 1/a/ 2)), (2n + l)a; <y< (2n + 2 )uj 

(A.33) 

where y = x — F is the hypergeometric function and 


= 


711+1/1 + 8 ^ I , P — n 



The coefficients 5+) are related to each other by 

= T 


7 = 3/4, n = 0,l,2,... 

(A.34) 


T = 


p^{n+l) 

\/2 cos7r(a — /3) 

87rr(7)^ 


5+) 


87rr(2-7)^ 


r(l-a)r(l-/3)r(l+a-7)r(l+/3-7) 
\/2cos7r(Q! — (3) 


(A.35) 

_A / V r>r\SD 'TTl r\, _ !•< \ I 

r(a)r(/3)r(7-a)r(7-/3) 

Using this recursion relation all the coefficients can be given in terms of 

The eigenvalues of the matrix T give the characteristic exponents as a function of 

. = gVX 


fi{s) = 


2^ 

r(i/ 4 )^ 


In < ^2 


cos 


tta/I + 8s 


± 


cos 


7rVl + 8s' 


(A.36) 
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The maximal value is obtained with values s = 2n^ — 1/8, where n = 1, 2,..and is 


f^max 


20F 

r(i/4)2 


ln(l + V2)^ 0.2377 


(A.37) 


The initial coefficients A®, are given in terms of initial conditions as 


\ - t( 

_V2_ _1_ 

T = -7r3/2 I r(7)r(i+o-7)r(i+/3-7) 2r(7)r(i-o)r(i-/3) 

I_ y/2 _1_ 

V r(2-7)r(a)r(/3) 2r(2-7)r(7-a)r(7-/3) 

The initial values are obtained from the values they have at the end of inflation. Any 
massive held during inhation behaves as 





r(z/.) 1 

r( 3 / 2 ) 


(A.39) 


where = 3/2 + e — t]^, where e and rj^ = are slow-roll parameters, 

is a Hankel function and the second expression applies to large scales. This results 
into = Tiiya — 1/2)Xct^ at large scales and therefore it also obeys the condition 
for getting rid of the non-locality. 
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